Theory of wave-front reversal of short pulses in dynamically-tuned zero-gap periodic 

systems 
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Recently we have shown that the wave- front of short pulses can be accurately and efficiently 
reversed by use of simple one-dimensional zero-gap photonic crystals. In this Article, we describe 
the analytical approach in detail, and discuss specific structures and modulation techniques as well 
as the required steps for achieving complete time-reversal. We also show that our scheme is only 
very weakly sensitive to material losses and dispersion. 
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PACS numbers: 42.25.Bs, 42.70.Qs, 42.79.Hp 



I. INTRODUCTION 

Time-reversal is one of the most spectacular yet elu- 
sive wave phenomena. A time-reversed pulse evolves as 
if time runs backward, thus eliminating any distortions 
or scattering that occurred at earlier times. This enables 
light detection, imaging and focusing through complex 
media (see e.g., [2-1^1 and references therein) with appli- 
cations in diverse fields such as medical ultrasound 0- 
0] , communication systems and adaptive optics [B43| : su- 
perlensing Q, ultrafast plasmonics Q, biological imag- 
ing [loj , THz imaging and quantum information and 
computing jl^ . 

For low frequency waves (e.g., in acoustics, spin waves, 
elastic waves etc.), time-reversal can be accomplished by 
electronic sampling, recording, and playing back [2|-|j]. 
This is possible since in this frequency range, the pulse os- 
cillates on a scale slower than electronic sampling speed. 
On the other hand, for high frequency waves, specifically, 
optical electromagnetic waves, the period of the pulse os- 
cillations (namely, the inverse of the carrier frequency) is 
too fast to be sampled properly by even the fastest elec- 
tronic detector. Moreover, for sufficiently short pulses, 
even the pulse's envelope is too short to be sampled elec- 
tronically. 

The standard solution in the optical regime is to 
use nonlinear processes such as Three- Wave Mixing or 
(nearly-degenerate) Four- Wave Mixing (TWM or FWM, 
respectively) [ij] ■ When these techniques are applied 
to monochromatic (CW) beams (i.e., when the electric 
field of the pulse is given hy E = A{x, y, z)e~'^^°^ -\- c.c), 
a process known as Phase-Conjugation (PC), the inci- 
dent photon interacts with two CW pump beams of the 
same frequency such that the output photon has a neg- 
ative frequency, namely, the conversion is wq — > — wo- 
Mathematically, this is equivalent to a conjugation of the 
envelope A{x, y, z), hence, the source of the terminology. 
The wave is then accurately time-reversed. 

Application of PC with CW pump beams to a pulse 
{E = A{x, y, z, t)e~^^°* -\- c.c.) yields the same frequency 
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conversion, wg — — wq. In this case, however, the 
technique, which is more appropriately referred to as 
Temporal PC (TPC) [Hi, does not lead to a perfect 
time-reversal. Indeed, since every frequency compo- 
nent is shifted by the same amount, the envelope con- 
jugation is accompanied by an inversion of the spec- 
tral envelope with respect to the carrier frequency (i.e., 
A{uj — Wo) — )■ A*{—uj -\- Wo), with A representing the 
Fourier-transform of the pulse envelope). This process 
thus involves a non-zero momentum changes, which, in 
turn, cause phase-mismatch problems leading to lower 
efficiency and accuracy; eventually, the ability to re- 
verse pulses with a broadband spectrum is limited. More 
simply put, the reversal achieved by TPC is imperfect 
since although the time envelope {A(x,y,z,t)) is con- 
jugated, it is not time-reversed {t not inverted to —t). 
As a consequence, while TPC compensates for effects 
such as group- velocity dispersion and Self-phase modula- 
tion (which are the dominant distortions experienced by 
pulses propagating in long-haul optical fibers), it does not 
correct for other effects (such as odd-order dispersion and 
self-steepening) [ij]; Those become important for suffi- 
ciently short (usually sub-picosecond) pulses or near a 
zero group-velocity-dispersion point. Thus, in practice, 
TPC can accurately reverse only sufficiently long (and 
weak) pulses for which the effects which are not compen- 
sated for are anyhow negligible. 

Among the few possible remedies to this problem, the 
most successful one is to employ Spectral Phase Conju- 
gation (SPC), which is basically a phase-conjugation of 
the Fourier-transform of the pulse (i.e., A{ijj) — > A*{ijj). 
SPC was demonstrated using two cascaded TWM pro- 
cesses, see e.g., [H, [iBl and references therein. However, 
the complexity of the scheme [l^ makes it unlikely to 
evolve into a commercial product. A more recent imple- 
mentation of SPC using two time-lens fiber systems and 
FWM in Silicon waveguides has been demonstrated for 
picosecond pulses propagating in optical fibers [TtI . 

Several years ago, two additional time-reversal schemes 
have been suggested which rely on dynamically-tuned 
coupled resonators optical waveguide (CROW) struc- 
tures. The first, by M. Yanik and S. Fan (see [Tsj). 
involves an adiabatic modulation of the index of refrac- 
tion of two anti-symmetric CROWs. The second scheme. 
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suggested by S. Longhi [l^, is based on a time-gated 
spatial modulation of the index of refraction of the res- 
onators. Both schemes can be implemented with linear 
modulators, e.g., using the electro-optical effect or car- 
rier injection/depletion. They can be applied to pulses 
which are spectrally narrower than the bandwidth of the 
CROW. This, typically, corresponds to pulses not shorter 
than several picoseconds. To date, both schemes were not 
tested experimentally. 

Importantly, we note that both schemes [l^ [3l Per- 
form envelope (or wave-front) reversal (i.e., A(t)e~^'^°* — >■ 
A{—t)e~^'^°*) rather than complete time-reversal, which 
requires also the conjugation of the temporal envelope. 
Following the discussion above, one can understand that 
a hybrid scheme, based on a combination of envelope 
reversal and TPC (which conjugates the temporal en- 
velope) will yield perfect reversal, in a similar manner 
to However, as noted, the CROW-bascd envelope 
reversal schemes are spectrally limited, so the associated 
hybrid scheme can still reverse only relatively long pulses. 

In a recent study we have introduced an alter- 
native wave-front reversal scheme that requires a con- 
ceptually simple structure, namely, a dynamically-tuned 
zero-gap photonic crystal (PhC). When applied to op- 
tical pulses which are not shorter than several picosec- 
onds, a 100% reversal efficiency can be easily reached 
using linear modulators. Implementation of the scheme 
for shorter pulses is also possible, however, it requires 
a nonlinear mechanism such as Cross-Phase-Modulation 
and has a reduced efhciency. In our approach, there are 
no special phase-matching considerations that have to be 
taken into account, it can be implemented on-chip, and 
the required energies are significantly lower than with 
techniques which are based on a series of nonlinear wave- 
mixing processes. Furthermore, the scheme is, in prin- 
ciple, not limited to one-dimensional (i.e., optical fiber) 
systems, and can be implemented essentially in any wave 
system where the wave velocity varies periodically and 
can be modulated in time. 

While in [ij the scheme was described intuitively and 
demonstrated numerically, the analysis was described 
only briefiy. This article is dedicated to the detailed 
description of the theoretical approach and the imple- 
mentation details such as choice of materials, modula- 
tion techniques and the roles of dispersion and absorp- 
tion. In Section |lT] we recall the principles of broadband 
wave-front reversal in zero-gap periodic systems. In Sec- 
tions mil and IIVI we derive the envelope equations and 
in Section |V] we show how to solve them analytically in 
the weak-coupling limit. Section I VI I provides a discus- 
sion on implementation considerations and Section IVIII 
describes the numerical results. Explicit formulae for the 
coefficients are derived in the Appendices. 



II. PRINCIPLES OF WAVE-FRONT REVERSAL 
IN ZERO-GAP SYSTEMS 

Wave propagation in periodic systems has been the 
focal point of countless studies in several branches of 
physics. Most prominent of these systems are crystalline 
solids. However, a variety of artificial analogues are 
also extremely well studied. Among those are as pho- 
tonic crystals (PhCs) [1^ in the context of electromag- 
netics, phononic crystals [2 l| in the context of acous- 
tics, magnonic crystals [22|, |23| in the context of spin 
waves etc.. The study of such systems was almost ex- 
clusively associated with the forbidden gaps, which are 
energy /frequency regimes where waves cannot propagate. 
In particular, these gaps are the key to understanding the 
electronic and optical properties of solid-state systems; 
their artificial analogues have a variety of possible appli- 
cations, ranging from low-threshold lasers, wave guiding, 
storing, filtering and switching j20j to medical ultrasound 
and nondestructive testing |2l| . 

Although the occurrence of bandgaps is generic to pe- 
riodic systems, in some special cases, the gaps can have 
a zero width. Somewhat similarly, some materials and 
structures support two bands that cross symmetrically 
without forming a gap (see e.g. Fig. [2]). One of the most 
notable of such structures is graphene [l^] or its photonic 
crystal analogue (20j which is a two-dimensional periodic 
system with hexagonal symmetry. Other examples are 
chiral metamaterials (25j , transmission- line systems [26| , 
bi-axial crystals [l^] and spin systems [2^. We refer to 
such systems with a perfectly symmetric crossing as self- 
complementary systems, since they support wave propa- 
gation in opposite directions with the same momentum. 
Earlier studies of self-complementary systems have al- 
ready shown some unusual properties [l^, [H, [13, [29j . 
However, these systems are still relatively unexplored. 

In this Article, we show how to exploit the unusual 
band-structure of a periodic self-complementary system, 
namely, a PhC with a zero-width gap, for the purpose 
of wave-front reversal of broadband pulses. As described 
above, since the main challenge is to reverse pulses of op- 
tical frequencies, in this Article, we focus on time-reversal 
of electromagnetic pulses. 

The zero-gap system has two major advantages over 
a "standard" finite-gap system. First, while a pulse in- 
cident on a finite-gap system will mostly be refiected, 
a pulse incident on a zero-gap system would be almost 
perfectly admitted, with only frequencies overlapping 
distant gaps being repelled. The second advantage of 
the zero-gap system is that in contrast to the standard 
approach of time-reversal, which, as described in Sec- 
tion U is usually realized with a frequency conversion 
process based on TWM/FWM, the self-complementary 
system allows for a different approach, namely, via dy- 
namic tuning of the refractive index. Indeed, dynamic 
tuning of PhCs has been one of the recent hot topics 
in PhC research, with various applications such as fre- 
quency shifting, light switching, slowing and stopping 
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as well as many others studied theoretically and demon- 
strated experimentally, see, e.g., [sol - fssj . In the current 
context, the dynamic tuning enables the vertical tran- 
sitions between the positive and negative group-velocity 
bands. The proximity of the crossing bands allows for 
efficient transfer of energy between the bands using weak 
and slow modulations of the wave velocity. In particu- 
lar, the modulation can now be much slower than the pe- 
riod of the pulse oscillations so that the wave- front rever- 
sal can be realized electronically, e.g., using the electro- 
optic effect or carrier injection/depletion. However, the 
modulation should still fulfill several conditions. First, it 
should be much faster than the pulse duration or equiv- 
alently, the spectral content of the modulation should be 
much wider than that of the pulse. In such a case, the 
spectral content of the pulse is effectively constant, so 
that all the frequencies in the pulse are converted with 
the same efficiency. Only such non-adiabatic modulation 
ensures a frequency- independent frequency conversion re- 
sulting in a reversed pulse which is an accurate replica of 
the forward-propagating pulse. In that regard, we note 
that although the modulation can have fast spectral com- 
ponents, the possibility of transitions to higher bands in 
a zero-gap system is small, thus, enabling non-adiabatic 
modulation without a deterioration in performance. A 
second requirement is that the modulation be periodic 
in order to avoid any wavevector mixing, resulting in a 
vertical frequency-conversion, see Fig. [2ljb). 

Following [l| , here we analyze a simple one-dimensional 
zero-gap systems and confirm the intuitive arguments 
above. However, we emphasize that the ideas and tech- 
niques we use here can be employed also in other zero- 
gap systems, in one, two or three dimensions, as well as 
in other self-complementary systems (23 - [28| . 



then, the Maxwell equations reduce to 

d,E{x,t) = ^^lo^tH{x,t), (2a) 
d,H{x,t) = -eodt[n''{x,t)Eix,t)], (2b) 



and eventually to 



d,^E{x,t) - —dtt [n^{x,t)E{x,t) 



(3) 



which is a one-dimensional wave equation with a time- 
dependent velocity. 



"l 


"2 







FIG. 1: (Color online) Geometry of the pulse propagation 
through a layered PhC. 

We proceed with an analysis similar to that first per- 
formed in [3^ [131 in the context of pulse propagation 
in PhCs with a Kerr response [s^]. We adopt the uni- 
directional field formulation by defining 

E{x,t) = n-i{x,t)[F{x,t) + B{x,t)] , (4a) 
ZoH{x,t) = n^x,t)[F{x,t)-B{x,t)], (4b) 

where F and B are scalar functions representing the for- 
ward and backward fluxes in the PhC [s^ [s^. Substi- 
tuting in Eqs. ^ gives 



III. DERIVATION OF THE UNI-DIRECTIONAL 
WAVE EQUATIONS 

Consider an electromagnetic plane-wave pulse nor- 
mally incident on a ID Photonic Crystal (PhC) (see 
Fig. [T]) which is time-modulated in the following man- 
ner 

n(.T, t) = nphc{x) + Map{x)m{t), (1) 

with max{m{t)) ~ 1 and max{p{x)) = 1. Here, both 
the static part of the refractive index, nphcix), and 
the spatial profile of the modulation p{x) have period 
d. Such a periodic modulation docs not introduce any 
wave- vector mixing, hence, allowing for an accurate re- 
versal (see e.g., [131 )• The Maxwell equations in this case 
can be written as 

VxE{x,t) = -HodtH, 
V X H{x,t) = eodt n^{x,t)E{x,t) 
For a linear polarization at normal incidence, 
E = E{x, t)z, H ^ H{x, t)y, 



„ , , n(x, t) ^ TLf „ 

F^{x,t) + ^^Ft + —F 
c c 

„ , , n(x, t) „ TLf ^ 

B,{x,t)-^^Bt-—B 
c c 



2 V c n / 

1 + I^) F (5b) 

2 V c n / 



Note that the definition of the fiuxes F and B is valid re- 
gardless of the specific form of the refractive index distri- 
bution n{x,t) and without assuming any specific modal 
(spatial) basis. Also note that since the wave is assumed 
to be normally incident on the PhC, the fields are essen- 
tially scalar and the Maxwell divergence equations are 
automatically satisfied. For oblique incidence, or for a 
higher-dimensional PhCs, the problem requires a gener- 
alized treatment, see 1371. 



IV. DERIVATION OF THE ENVELOPE 
EQUATIONS 



We now define 



(6) 
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so that the system ([5]) can be written in matrix form as 
RWt=MK, (J) 

where 



where 



M 



fnphcix) 
V nphc{x)^ 

-f(l^ + lf) 



Mom{t) 



p{x) 




nt 



(8b) 



The static index variations, nphc, and the dynamical 
variations, p, vary spatially on a scale comparable to 
A; correspondingly, the fluxes vary in space on a similar 
scale. In addition, they vary in time on the scale of the 
period T. On the other hand, the pulse's envelope and 
the modulation m(<) vary on slower scales, Tp and Tmod, 
respectively. Thus, the solution ([SJ varies on more than 
one temporal and spatial scale. Accordingly, we adopt a 
multiple-scales approach (39| by assuming that the varia- 
tions on the fastest scale, defined as xq, to, correspond to 
variables independent from the variables describing the 
variations on the slower scales. However, in contrast to 
the standard multiple scales approach, where there are 
only two time scales (e.g., the period T and the pulse 
duration Tp] see [s^ or more specifically, [1^), the exter- 
nal modulation introduces a third independent time scale, 
Tjnod- In order to avoid this complication to the analysis, 
let us lump all the slow changes into one spatial and one 
temporal scale, defined as xi and ^i, respectively. Then, 
Eq. ([T]) becomes 

i Rp^c^xo) + MoR^^{xo)m{ti), (9) 

and we can make the substitution 

dt df„ + df, , d.j; d^^„ + d^, , (10) 

so that the matrix M becomes 

= \ -^dt^n- ^dxoU c{dxo + dxi) - dt^n 

(11) 

In order to complete the scale separation, let us rewrite 
the terms involving the spatial derivatives in M (|TT|l as 
a sum of a static and dynamic parts. Using Eq. ©, it 
can be shown that 



n 



'PhC 



(12) 



^ Mom{ti) {nphcp'jxg) - n'p,^^p{xo)) 
nphcixo) [nphcixo) + Mop{xo)m{ti)] ' 

where the sign ' stands for a derivative with respect to the 
argument, a notation used only when it is unambiguous. 

We can now separate M into the different order con- 
tributions as 



(13) 



/ 1 "Phc \ 



2 lie 



(8a) M}'\xi,ti;xo) = 



c ( /3(xo,ti) 



P{xqM) 
Vdx, - MoC^m\t,)p{xo), 



14b) 



and where we have defined the space and time-dependent 
function 



l3{xo,ti) = m{ti)- 



nphcixo)p'{xo) ~ n'p^^p{xo) 



nphc{xo) {nphc{xo) + Mop{xo)m{ti)] ' 

(15) 

and the constant matrices 



V 



c 
-cl ' 



and 



The matrix M*-*^-* is the operator of the unperturbed sys- 
tem (i.e., in the absence of time- modulation) while the 
matrix m}^\ together with the modulation term on the 
left-hand-side of Eq. (or equivalently, the second term 
in Eq. (HI), describe the slower dynamics, i.e., the evo- 
lution of the pulse envelope (wave-front) and the effects 
of the modulation; the dependence of M*-^-* on the fast 
scale xq is removed by integration, see below. 

The formulation presented so far is valid for any ID 
periodic optical systems (or PhCs). In what follows, we 
limit the discussion to periodic systems in which two 
bands have minimal interaction, leading to a gap with 
effective zero width, or effectively, a crossing of the two 
bands. Such crossing points typically appear between 
the higher bands of PhC structures. In most cases, they 
are not perfectly symmetrical with respect to the cross- 
ing point, however, nearly- or even perfectly-symmetric 
crossings can be found. Probably the simplest example 
of a perfectly symmetric zero-gap periodic system is a ID 
layered PhC for which the indices and thicknesses of the 
two layers satisfy 



nidi = ^2^2/5, 



s = 1,2,3,. 



(16) 



In such systems, the s + 1 gap at vacuum wavelength [5£ 



Ac ~ 2nidi = 2n2(i2/s, 



(17) 



has exactly zero width. For dispersive materials, i.e., 
when the refractive indices depend on the wavelength, 
the QWS condition ([TB)) has to be satisfied exactly at 
the crossing point in order for the gap to maintain a 
zero-width. For simplicity, we assume that the disper- 
sion is negligible, and defer the discussion of the role of 
dispersion to Section fVl CI 
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(b) 




V < 



° K(;t/d) '' K^ = K 

FIG. 2: (Color online) (a) Dispersion relations of a QWS. (b) 
Schematic diagram of the reversal in the frequency domain. 
Only the part of the dispersion curve close to the crossing 
point Kc = is shown. 

When s = 1 , this system is called a quarter- wave stack 
(QWS) PhC (see e.g., liO,!!^). For simphcity, in what 
follows, we analyze the case of a zero-gap in the QWS 
system, sec Fig. [^^b). As shown in Appendix iBl PhCs 
with other values of s yield the same coefficients, hence, 
the same reversal efficiencies. 

Near a gap. the solution of Eq. ([7]) consists predom- 
inantly of two spectral components, one on each of the 
crossing bands. Each can be written as the product of 
a carrier uni-directional Floquet-Bloch (FB) mode ^ f /b 
on either the positive (/) or negative (&) group-velocity 
band with a slowly varying envelope i.e., 

W ^ [/(xi,ti)*^(xo) + 6(a.i,ti)i^fc(a;o)]e-'"=*°+c.c., 

(18) 

where ujc = 27rc/Ac is the band crossing frequency 
and c.c. stands for the complex conjugate. The uni- 
directional FB modes are the eigensolutions of 
Eq. (O, i.e., they satisfy 

- iuJanphc{xQ)%^{xo) = M'°n-^0)*a, (W) 



where = Wj/f, is the central frequency of the pulses on 
each band, chosen at the same FB wavevcctor. 

In what follows, we neglect the complex conjugate 
terms in Eq. (|T8|) . This is justified for pulses which are 
sufficiently narrow spectrally, i.e., for pulses for which the 
separation of scales (|10p is meaningful. We also ignore 
the population of any other distant band [11] . Indeed, as 
shown in [s^, the population of distant bands (or com- 
panion components) is manifested by the appearance of 
a non-zero second order term in the k ■ p expansion of 
the band structure. This term accounts for the group- 
velocity dispersion (GVD) whose magnitude scales in- 
versely with the temporal duration of the pulse. Since, 
however, the QWS system exhibits strictly linear disper- 
sion in the vicinity of the band crossing, GVD vanishes, 
so that population of distant bands does not play a role 
in the dynamics even for very short pulses [6 Of . 

The central frequency of the forward propagating 
input pulse can be detuned from the crossing point 
{Kc = 0,u!c)- Since the band structure is symmetric 
with respect to the crossing point (see Appendix |X] and 
Fig. [2]), and since the frequency conversion is vertical (see 
Fig.djb) and discussion after Eq. we can define 



Sujf — LOc — ujf = —SuJb = (20) 

We also assume that the detuning is small compared with 
the wave period, A <C 2tt/T, in agreement with the above 
assumptions on the separation of scales (fTO| . 

Substituting Eq. dTB]) in Eq. (O, muhiplying by *J 
and respectively, and integrating over xq allows us 
to remove the dependence on the fast scale xq. In par- 
ticular, this leads to the following equations for the time 
evolution of the forward/backward envelope components 



(l + Momiti)mf}^ ft, ixi,h) + Vfff,, - (zWcMo™//m(ti) - Mom'mf} + iSuf + Mo(3ff{h)'^ f (21a) 

(l + AfoTO(ii)mi°^) btAxi,ti) + Vbbb^, - (iu;,Mom["^m(ii) - Mom'mH^ + iSuJb + MoMh)^ b (21b) 

= -Momf^m{ti)ft, - Vbffxi + ^^o {iiOcmf^ m{ti) - m'mi^ + Pbfiti)^ /, 



I 

where we defined the coefficients as the integrals over the uni-directional FB modes, as follows. 



Vab / :^Uxo)V%ixo)dxo, (22a) 







(0) 

Kb 



ab 



/ ■^i{xo)Rjxo)Mxo)dx^, (22b) 
Jo 



d 



*I(-^o)n (xo)C^*fc(xo)dxo, (22c) 







Pab{ti) = - / ±i{xo)p{xQ,ti) ^b{xo)dxo. {22A) 
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Eqs. (|21a[) - ()21b|) arc transport equations with time- 
dependent velocities, coupled by the modulation m(ii). 
The equations for / and b are symmetrical, however, the 
initial condition is highly asymmetric, consisting exclu- 
sively of forward waves (see also Section IVlip . Moreover, 
note that the form of the coefficient of the dynamic terms 
/tj and bt^ (namely, I + Mq...) implies that a modulation 
that increases the refractive index {Mq > 0) will not 
give the same results as a modulation that decreases the 
refractive index {Mq < 0). This effect, however, is signif- 
icant only once the change in the pulse velocity due to 
the modulation is appreciable. 

The existence of simple explicit analytical formulas for 
the FB modes of the layered system (see e.g., [lO]), allows 
us to compute the coefficients (|22ap - (|22cp analytically. 
Near the band crossing, the obtained expressions are rel- 
atively simple, and shed light on the physical meaning of 
the various terms in Eqs. (|21ap - (j21b[) . First, the diago- 
nal velocity coefficients Vfj (vhb) represent the velocit y of 
the forward (backward) pulses. Indeed, it is shown in [36[ 
that the Vaa term is equal to the first term in the k ■ p ex- 
pansion of the band structure, i.e., to the group- velocity; 
also, in Appendix IB H it is shown that 



-Vbb 



duj 
'dk 



(23) 



The off-diagonal velocity coefficients Vat are responsi- 
ble for reflection coupling (i.e., for horizontal transitions 
k ^ —k) |3f|. Intuitively, the strictly linear dispersion 
relations at the crossing point imply that no reflection 
coupling is taking place, i.e., that Vh f = Vfb = 0- Indeed, 
this is confirmed in Appendix IB II |61| . The vanishing 
of the reflection coupling ensures that the modulation- 
induced transitions between the bands are strictly verti- 
cal so that the frequency conversion is free of wavevector 
mixing, see Fig.^b). 

The diagonal mia^ coefficients describe the phase accu- 
mulation induced by the modulation. Indeed, it is shown 
in Appcndix lB 2l that these diagonal coefficients are iden- 



tical, i.e. 



md 



- (0) 

^ff 



'fob 



and that for a piecewise- 



uniform modulation, they are given by 



1 f n.mi ^ n„i2 
ni 712 



(24) 



where rimj is the amount by which layer j is modu- 
lated [Hi. 



The diagonal rn^a coefficients represent the change of 
the electromagnetic energy induced by the modulation. 
In Appendix IB 31 it is shown that for a piecewise- uniform 



modulation, the coefficients m^aa and m^aa are equal. 



The off-diagonal coefficients = m^y are respon- 
sible for reversal coupling, i.e., to the vertical transitions 
uj + Soj ^ u! — 6uj, see Fig. EJb). One can interpret these 
coefficients as the measure of the imperfect spatial over- 
lap (or cquivalently, the phase mismatch) between the 
modes immediately above and below the crossing fre- 
quency. Their mathematical form (sec Eqs. (j22bp - (j22c|) ). 
which is reminiscent of the orthogonality integral (see 
e.g., Eq. (|A1[) ). implies that their size should increase 
with the index contrast between the layers. Indeed, it is 
shown in Appcndiccs IB 2| | B 3l that for a piecewise-uniform 
modulation and zero-detuning, the coupling coefficients 
arc given by 



Thus, the reversal coupling tends to zero as the index 
contrast approaches zero or if the modulation maintains 
the QWS condition p^ . i.e., if no gap is opened due to 
the modulation. In all other cases, the dependence of 
the coupling coefficient on the indices is non-trivial, see 
also Section fVl Bl It can also be shown that for non-zero 
detuning, mod becomes complex. 

Lastly, the /3ab coefficients (j22d| are particularly in- 
triguing. They arc defined in terms of an integration 
over /? (|15p . For the layered structure under consider- 
ation, n'pf^fj (and as well for a piecewise continuous 
modulation) vanishes everywhere except at the interfaces 
between the layers, where it becomes inflnite. In these 
cases, the coefficients (|22dp can be evaluated by smooth- 
ing the index discontinuities over a finite length L, com- 
puting the FB modes numerically, and then taking the 
limit i — >■ 0. Our numerical computations, performed for 
several configurations of indices, thicknesses, and sizes of 
time-localized modulations, indicate that for a layered 
PhC, these terms vanish, jiab{L — )■ 0) = 0. Neverthe- 
less, any other configuration (namely, a smooth index 
PhC, spatially non-uniform modulations etc.) requires 
the speciflc numerical evaluation of j3ab- 

Due to all the above, for a piecewise-uniform modula- 
tion, Eqs. ((2Ta|) - (j21bp simplify into 



(1 + MomdTO(ti)) /ti(a-i,ti) + Vgfx^ + Momdim' ~ iujcm{ti)) f ~ iSujff 

= Mom'-°J {iuj,m{t,)b~m'b-m{ti)btj, (26a) 

{I + Momdm{ti))btj^{xi,ti) - Vgb^^ + Momd{m' - iujcm{ti))b - iSujbb 

= MomlaiiujMti)/ - m' f - m{ti)ft,) . (26b) 
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In what follows, wc focus on this set of equations. 

Before finishing this Section, we would like to comment 
on alternative derivations of envelope equations for pulse 
(or beam) propagation in periodic structures. First, simi- 
lar equations to Eqs. (|26aP - (j26b[) can be derived using the 
Slowly- Varying Envelope Approximation (SVEA), i.e., 
by assuming that the carrier wave is a Fourier mode (i.e., 
a plane-wave) rather than a Floquct-Bloch (FB) mode. 
However, not only the uni-directional wave formulation is 
simpler, more elegant and more accurate than a deriva- 
tion base on the SVEA [4l|, but for PhCs, the SVEA 



WEAK COUPLING ANALYSIS 



is not strictly valid, since in the presence of A-fast index 
variations, the second derivatives are not small compared 
to the other terms. An exception is, possibly, the case 
in which the index contrast is small (or in other words, 
that the grating is shallow), sec e.g. [38|. Indeed, in that 
case, the difference between the Fourier and FB bases is 
small. 

A different approach described by Craster et al. is to 
derive envelope equations for the wave equation (rather 
than for the directional flux equations, as in our deriva- 
tion). This approach was employed to a large variety 
of periodic structures in one and two dimensions, see 
e.g., [H. 

Finally, a different set of carrier modes, computed from 
the transmission and reflection data of a given finite size 
PhCs, were used in coupled Helmholtz equations in the 
context of three quasi-monochromatic wave mixing in 
PhCs, see [i^ and references therein. However, this ap- 
proach is only suitable for pulses much longer than the 
PhC, since in that case, the reflections from the ends of 
the PhC system are important. 



In principle, Eqs. p6a|) - (|26bp are significantly simpler 
to solve numerically compared with the wave equation ([3]) 
since they do not require to resolve the A— and T— scale 
oscillations. In practice, however, the numerical solution 
of the envelope equations is non-trivial (see Section IVIII 
for details); Moreover, they are also difficult to solve an- 
alytically. In order to facilitate the analysis, recall that 
in our case, only the forward component (i.e., the posi- 
tive group-velocity band) is initially populated. Thus, as 
long as the backward component is small compared with 
the forward component, we can neglect all the coupling 
terms in the equation for the forward wave (|26ap . i.e., we 
can set its right-hand-side to zero [g^. In other words, 
we assume that there is only forward-to-backward cou- 
pling, or simply, that the coupling is weak. In this case, a 
closed form solution can be obtained to Eqs. (|26ap - (j26b|) . 

In order to derive that solution, note that the charac- 
teristics of Eqs. (|26ap - (|26bp arc solutions of 



dxi = ± 
We integrate and define 



Vqdtl 



1 + Momdm{ti) ' 



tl 



1 + Momd'm{t[) 



This is equivalent to transforming into frames moving 
with each of the pulses. Then, Eqs. (|26ap - (|26bp reduce 
to 



{1 + Momamih)) ft, (^xY\h^ 



{iuJcMQ'm{ti)'md + iScof — M^nidm') /, (27) 

{iuJcMom{ti)'md + iSuJi, — M^mdm') b (28) 

/ Z"*^ dt' ^ 

Mom(ti)mod (/tl - i(^cf) - Mom'modf i xf^ - "ivg / , , — j-tt, h 

V J-oo 1 + Momdm(t-^) ^ 



We now define t = (Mom^Wc) 
Eq. (P7)) can be written as 



The solution of 



f{x^/\h) ^ ef'-^^'^*'^^<fo{x^/^) 



^(t'l) 



+ iSujf — M^mdm' {t'l) 



1 



, (29) 



MQmdm{t\) 

where /o is the envelope of the pulse traveling in the for- 
ward direction before the onset of the modulation. Sim- 
ilarly, we define 



l + Momdm[t\) 



(30) 



r 



Note that for zero detuning and uniform modulation, 
= Pf- Also note that the definitions of and 
Pb confirm the physical interpretation of the coeffi- 
cients ([24]), namely, the phase accumulation and change 
of electromagnetic energy due to the modulation, as de- 
scribed in Section HYI Hence, 
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*i mit'Aft' -iuJcm(t\)f + m'it\)f X , 



1 + Momdm{t[) 

We substitute Eqs. ((27)) and (|29l) into Eq. ([3T|) and get after some algebra 



(4' 



1 + Momdm{t'{) 

(32) 



where 



h{t[;Mo,ujf) 



[I + Momdm{t[)f 



(33) 



At times long after the modulation has ended we can set 
the upper limit of both integrations to oo. Thus, Eq. ([5^ 
shows that the backward wave is given by a convolution 
of the forward wave with h(ti] Mq, ujf), which has the role 
of the impulse response of the system. This mathematical 
form of the backward wave is reminiscent to the results 
of wave-mixing-based spectral phase conjugation [T5| in 
which the backward wave is a given by a convolution of 
the signal with both pump waves. 

It follows from Eqs. p2p -p3 p that accurate reversal can 
be achieved only in the non-adiabatic limit, i.e., when 

Trnod < Tp, (34) 

in which case the impulse function is essentially a delta 



function. 

As an example, let us assume a Gaussian modulation 



m(t) = e "^mod 



(35) 



For such a symmetric modulation, one expects the contri- 
bution of the m'{t) term in h (j33]) to be small compared 
with the contribution from the ujfmit) term. Indeed, 
first, m'[t) ^ T^^ld ^ '^f- Second, m'{t) has odd symme- 
try, so that contributions from times before and after the 
middle of the modulation to nearly cancel. This was eon- 
firmed with asymptotic calculations (not shown). Thus, 
in what follows, we neglect the m'{t) term in h ([55)) . 

For a small modulation, we can also neglect the Mq 
term from the denominators, i.e., we effectively assume 
that the pulses' velocities are not affected by the modu- 
lation. In this case, x[^'''^ = xi =F Vgti. Using Eq. (|30)) . it 
then can be shown that Eq. (|5^ reduces to 



b(xf\t,) = i^;Afom„de''"''*^-*^«'"''™(*'i)+^/-~'"(*f H m{t[)e'^''^ fo (x'{'^ ~ 2vgt[) dt[. 

^ ^ J — OO 

(36) 



This integral can be solved exactly for a Gaussian pulse 



fo 



VgTp 



e s p , 



(37) 



where Tp = 0.7W/vg, with W being the full-width at 
half-maximum (FWHM) of the incident pulse. In this 
case, the convolution integral is 



with 



eff mod p 



'dt[, 

(38) 
(39) 



Then, it can be shown that 



VT^Tefffo 



with $ = 



T^ + '^Tl^d^ ). This shows 



that the envelope of the backward pulse is indeed a prop- 
erly reversed version of the envelope of the forward pulse; 
however, it is also broader than the forward pulse due to 
the convolution with the modulation [g^. More gener- 
ally, arbitrary asymmetric pulses will also undergo some 
distortion due to that convolution. These effects can be 
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minimized in the non-adiabatic hmit (j34[) . where we get 



VgTp 



(40) 



Finahy, substituting in Eq. p6p gives an overah reversal 
amplitude of 



maxti\b{xi,ti)\ = ^/TrTrnodUJfMo\mod\e' 



7o(0). 



(41) 



Note that the solution is symmetric with respect to the 
sign of Mq, because we neglected the effect of the modu- 
lation on the pulse velocity. It also shows an exponential 
(Gaussian) decrease of the efficiency with the detuning. 

The calculation above can be generalized for 
arbitrarily-shaped pulses. Moreover, the calculation can 
be performed without the simplifications applied to the 
impulse response p3p . Inclusion of these terms yields 
0{Mq) corrections to the amplitude of the reversed pulse. 
However, such a calculation docs not take into account 
even 0{Mo) corrections to the amplitude and profile 
of the forward wave. This can be amended by self- 
consistency iterations, but this would be tedious and 
would probably extend the validity of the perturbation 
analysis only slightly. More generally, once the backward 
wave becomes comparable in magnitude to the forward 
wave, in particular, for sufficiently long modulation times 
and/or strong coupling, the weak coupling solution (|4ip 
is not valid anymore. In all those cases one needs to 
solve the envelope equations (|26ap - (|26bp without approx- 
imation. Nevertheless, despite these limitations, we show 
in Section IVIII that the weak coupling solution (|4T|) pro- 
vides a very good approximation to the actual dynamics 
(obtained by solving the wave Eq. ([3]) and the envelope 
equations (|26a[) - (|26b[) ) even for conversion efficiencies as 
high as 50%. 



VI. DESIGN RULES 

We now turn to discuss the practical aspects of design- 
ing our Reversal Mirror (RM) for electromagnetic pulses. 
We take into account realistic material parameters, dis- 
cuss efficiencies and a variety of other practical issues. 



A. RM thickness 

The length of the RM, D, should be chosen such that it 
contains all the pulse during the modulation (in a similar 
manner to 0, [31). As a simple estimate, this requires 
D > NW where W w VgTp/0.7 is the FWHM or equiva- 
lently. 



with typically, N > 3/0.7 ^ 5. Thus, for example, a few 
cm-long RM can contain and reverse pulses not longer 
than a few tens of picoseconds. 

The thicknesses of the layers determine the crossing 
frequency, hence, the carrier frequency, according to 
Eq. (|T7l) . Note, however, that the parameter s in Eq. (fT6)) 
allows to scale the thickness of one of the layers. 



B. Materials 

1. Modulation technique 

We distinguish between two classes of modulations. 
First, picoseconds and longer modulation times can be 
done externally, e.g., by using the electro-optical (Pock- 
els) effect |13j or using carrier injection/depletion in semi- 
conductor materials, see e.g., [33 • These techniques yield 
a modulation which is, to leading order, uniform; both 
can be realized with essentially linear modulators (see 
also [illillll). 

Large electro-optical coefficients arc found in ferroelec- 
tric materials such as LiNbO^, BaTiO^ or BaStTiO^, all 
which can be modulated on the GHz scale and faster. In- 
dex modulation via carrier injection/depletion can be ob- 
tained in a variety of semiconductors. For example, index 
modulations of the order of 1% at about lOOGHz were 
achieved in silicon devices [H, H^l . Even faster modula- 
tions of Ips can be achieved in ion-implanted silicon [45j . 
Gomparable modulations may be attained in other semi- 
conductors such as GaAs. 

Shorter modulation times can only be achieved using 
intense ultrashort pump pulses and rely on a nonlinear 
mechanism such as Gross-Phase modulation (XPM) or 
carrier injection/depletion [35ll46j through a pump-probe 
type experiment. In such cases, the modulations may not 
be uniform or even periodic, unless the geometry of the 
modulation is designed appropriately, see e.g., |47| . 

We note that in the case of XPM-based modulation, 
the RM structure need not only to have a periodically 
varying refractive index, but also a periodically varying 
cubic nonlinear (Kerr) coefficient. A systematic study of 
pulse propagation in structures with a periodically vary- 
ing (positive) Kerr coefficient was conducted by Sivan 
et al. [H, |4^ . This work was elaborated and extended 
to systems in which the magnitude and sign of the Kerr 
coefficient vary periodically, as well as to structures in 
which both refractive index and Kerr coefficient vary pe- 
riodically. For a recent review of light propagation in 
such structures, see review by Malomcd et al. [50|. 



Tp < Q.ID/Nvg 



2. Refractive index 

As shown in the formula for the coupling coeffi- 
cient (pS)) . high reversal efficiency can be achieved if (a) 
there is a high contrast between the refractive indices of 
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both layers and (b) the modulation is as different as pos- 
sible from such that maintains the QWS condition. Both 
conditions ensure that a gap as wide as possible is opened 
by the modulation. In particular, this implies that it is 
far better to modulate a single layer (e.g., by choosing 
one of the layers to be air), rather than two layers by 
the same amount. In this case, the coupling coefficient is 
given by 



\mod\ 



\ni - n2\ 
2nj{ni + ^2) ' 



(42) 



where j is the index of the modulated layer. In the limit 
of a high contrast, mod 1/2?^^ so that it is clearly 
advantageous to modulate the low index layer. If, on 
the other hand, one modulates the high index material, 
the refractive index of the high-index layer required for 
optimal reversal is ni^opt = ^2(1 + still, since in this 
case the dependence of the reversal coefficient on ni is 
rather weak, there is only a small reduction in efficiency 
for indices slightly different from rii^opt- 

These configurations can be realized with material 
combinations such as silicon or LiNbO^ and air, both 
giving rise to mi,f ^ 0.08. An alternative is to use a 
low- index ferroelectric with a high- index material (e.g., a 
semi-conductor) . Then, if only the ferroelectric is modu- 
lated, the coupling coefficient can have comparable mag- 
nitude. Another possibility is to use high refractive- 
index materials, either for constructing mid-IR (5l| and 
THz RMs. The associated efficiencies may be slightly 
higher. 

Finally, an even higher reversal efficiency can be at- 
tained if the refractive index of one layer is increased 
while the refractive index of the other is decreased. This 
is possible, for example, with XPM-based modulation in 
structures in which the Kerr coefficient has different sign 
in each layer (see [H^ and Section IVIB 1[) . 



the band structure calculation (see [40] or more specifi- 
cally. Appendix |^ . In this case, dispersion will cause a 
deviation from the strictly linear dispersion relations and 
will be manifested via the addition of high-order time- 
derivative terms (GVD and higher order dispersion) in 
each envelope equation. These will lead to broadening 
and even distortion (in the presence of significant higher- 
order dispersion) of the pulses, possibly to a different 
extent on either band. However, since the structural dis- 
persion of the PhC, dominates over material dispersion 
(which is a weak effect in dielectrics), dispersion is not 
expected to be significant, especially due to the short 
lengths of the proposed RMs. 



D. Losses 



So far, we have also ignored the possibility of losses in 
the RM (i.e., absorbance or scattering). As a rough esti- 
mate, if one of the layers has non-zero absorbance, then, 
it is sufficient for the absorption coefficient of the lossy 
layer, a = 2fcon", to satisfy a^D <C 1. For example, 
typical losses in silicon nanostructures (e.g., channel or 
ridge waveguides where losses originate mainly from scat- 
tering from the roughness of the interfaces) vary between 
0.05 - 1/cm m. For a silicon-air QWS (nst ~ 3.5, so 
that d ^ (ni -\- 712)^1/71.2 A.bdi), this corresponds to 
the condition a <C -jfjj ^ Thus, such losses can be 
sufficiently low even for a few cm-long RM. If the ma- 
terial of choice has even higher losses, one can choose a 
different value of s in order to reduce the thickness of the 
absorbing layer. Thus, overall, losses are not expected 
to cause any significant reduction in performance of our 
RM. 



C. Material dispersion 

In the derivation of the envelope equations we have 
ignored the effects of dispersion. These arc unavoidable, 
especially since our design prefers high index materials 
and long samples. Dispersion affects several aspects of 
the RM. First, since now the refractive index depends on 
the wavelength, it is clear that the QWS condition 
cannot be satisfied for all wavelengths. However, it can 
be easily shown using the analytical expressions for the 
band structure that in order to maintain a zero-gap, it 
is sufficient for the QWS condition (|16p to be satisfied 
exactly at the crossing wavelength ((TT)) . 

Second, dispersion will affect the dynamics of the 
pulses inside the RM. Self-consistent inclusion of dis- 
persion in Maxwell equations requires to compute the 
dynamics of the Polarization so that the refractive in- 
dex arises naturally from the constitutive relations [53| . 
However, for the purpose of the derivation of envelope 
equations, it is sufficient to introduce the dispersion into 



VII. NUMERICS 

In order to validate the analysis, we have performed 
extensive numerical simulations. First, the wave equa- 
tion ([3]) was solved with standard absorbing boundary 
conditions at both ends of the domain, namely, 

„ / V n(x = O.t) , 

E,{x = 0,t) = ^ U ~ 2E,nc.t) , 

c 

„ , „ X n{x = D,t) ^ 

c 

Second, the envelope equations (|26a[) - (|26b[) were solved, 
again, with standard absorbing boundary conditions for 
both envelopes with only the forward envelope initially 



excited, namely, 

/j;(.T = 0, = (/t - 2/i„c,t) , 

c 

c 

c 

Although the numerical solution of these equations 
does not require to resolve the T- (and A-) scale oscil- 
lations, it still poses a non-trivial challenge. First, note 
that since the modulation is uniform in space, then, un- 
like the case of light propagation through nonlinear Kerr 
media [s^ H^l where the modulation is caused by the 
moving pulse itself, Eqs. (|26ap - (j26bp cannot be reduced 
to a set of two coupled ordinary differential equations. 
For this reason, simp le numerical schemes such as the 
one described in [S^l are not applicable to the current 
problem and one must employ partial differential equa- 
tions techniques. Second, the absence of diffusion from 
the transport equations (|26a|) - (|26bp makes any discretiza- 
tion scheme unstable. In order to retain numerical stabil- 
ity, one can add small artificial diffusion, or alternatively, 
use the Galerkin-Ritz technique which is relatively non- 
diffusive. 



(a) (b) 




t(ps) 0.13 x(nm) 200 



FIG. 3: (Color online) (a) Amplitude of an asymmetric pulse 
(Eq. ([31; red line) and the associated backward wave envelope 
(Eq. (|26bp : blue line) at the input side of the QWS as a func- 
tion of time, (b) A spatio-temporal contour map of forward 
and backward envelopes in (a). 

In Fig. [31Ja), we plot the pulse profile at the input side 
of the PhC as a function of time. The reversed pulse has 
a somewhat lower-amplitude, but the leading and trail- 
ing edges have clearly exchanged roles, so that wave- front 
has been accurately reversed. We also show that the so- 
lutions of the wave equation ([3]) and the envelope equa- 
tions (|26a[) - (|26bp are in excellent agreement. Fig. EKb) 
shows a spatio-temporal contour map of the pulse prop- 
agation. It shows that the modulation, occurring once 
the pulse is in the middle of the RM, causes a fast and 
complicated dynamics after which the pulse splits into a 
reversed and a (somewhat delayed) forward component. 
Note that the modulation causes a delay of the forward 
propagating pulse (occurring even for very small modu- 
lations). This effect is irrelevant to the performance of 
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(a) (b) 




-5 5 -8 8 

10^ A 10^ A 



FIG. 4: (Color online) (a) Reversal of a Tp = lOOfs, unit- 
amplitude Gaussian input pulse (|37|1 under a Gaussian mod- 
ulation p5|) as a function of the index change A for ni — 3.45, 
712 = 1, Ac = 1550nm, and Tmod ~ lOfs. Shown are numerical 
solutions of the wave equation (Eq. ((3]); blue dots) vs. the so- 
lution of the envelope equations (Eq. (|26b|) : black circles) and 
the analytical solution (Eq. (|41|l : red solid line), (b) Same as 
(a) for Tmod = Ips and Tp = lOps. 



the RM, and is not captured by the analysis performed 
in this article. 

A comparison of the reversal efficiencies as a function 
of the modulation strength, which also includes the an- 
alytical approximation (|4ip , is shown in Fig. 2] We em- 
ploy realistic parameters corresponding to a Silicon-air 
QWS with Ac = 1550nm; such a PhC can be fabricated 
on-chip with current technology [35| . First, we perform 
simulations of ^ 20-cycle pulses. Fig. IDJa) verifies that 
all three solutions are in good agreement, thus, validat- 
ing the analysis. The reversal efficiencies in this case are 
rather low, however, they are comparable to efficiencies 
of nonlinear conversion processes of such short pulses. A 
possible solution is simply to amplify the signal after the 
reversal process. 

In Fig. [UJb) wc show simulations for longer pulses for 
which memory and running time required for the solu- 
tion of Eq. ([3]) are beyond standard available computing 
resources. Thus, we only show the solutions of the en- 
velope equations (|26ap - (|26b[) vs. the analytical solution. 
As in Fig. Hl^a), there is very good agreement between 
the numerical and analytical solutions up to high effi- 
ciencies (^ 50%). At even higher efficiencies, the analyt- 
ical approximation overestimates the reversal efficiency 
as given by the solution of the envelope equations. This 
is because at such high efficiencies, the forward wave am- 
plitude is significantly decreased, so that the forward- 
to-backward wave coupling decreases as well. Neverthe- 
less, the solutions of the envelope equations show that 
a 100% reversal efficiency can be achieved in this con- 
figuration with index changes only slightly higher than 
those predicted analytically. Thus, our RM has compara- 
ble performance to the previously suggested time-reversal 
schemes in index-modulated coupled-resonator arrays of 
optical waveguides [ll,[l3- However, in contrast to these 
schemes, the QWS admits pulses of very broad spectrum. 
In addition, our system does not suffer from a deteriora- 
tion of performance due to non-adiabatic modulations. 
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Accordingly, to 0{5lu'^) accuracy, for any value of s we 
get 



A = l + i doj — 

c 2nin2 

so that by Eq. (|X2|) . 



nidi (n2+ni)^ 



2nin2 



Appendix A: Floquet-Bloch modes near the crossing 

point 

The Floquct-Bloch modes are solutions of the 
Hclmholtz equation, obtained by Fourier transforming 
the wave equation ([3]) with Mq = 0. In a single unit cell, 
they are given by 



cos {K{uj)d) = 1 — 



nid 



-Suj 



(?i2 + ni^ 
2ni7i2 



Since we are interested in a transition between modes at 
a given (say positive) k, it follows that 



K(OJ) = — \ooj\ > 0, 



(A5) 



1 



■e"'^u(a;). 



" ^ ^ coe''=2(^-'^i) + doe"*''"(='~''i', di <x <d, 

where kj = ujnj{uj)/c, is a normalization constant 
chosen such that 



4>l{x)n^{x)(t)h{x)dx ^ 5ab, (Al) 



and the FB wavenumber is given by 

n[uj) = ^acos {Re [A]) , (A2) 

with A being the first element of the transfer matrix of 
the layered system, see [iOl and Eq. (|A3cP below. Note 
that for simplicity of notation, wc omit the k subscript 
from u, and ^P. 

For normal incidence, the coefficients for the first layer 
are given by [S, Eq. (6.2-7)] 



Oq 



B 



2 7iin2 



ha = e-""^-A, 



A 



-iki di 



COs(fc2C?2) 



I Tin + n 



2 ni7i2 



(A3a) 
(A3b) 
isin(fc2rf2) ;A3c) 



The coefficients for the second layer are given by [40, 
Eq. (6.1-11)] 



where the absolute value over the frequency ensures that 
the modes under consideration have the same Bloch mo- 
mentum, as required for a vertical transition, and the 
absolute value over the indices ensures that the group- 
velocity is real even for negative index materials. For 
simplicity, in what follows we drop the absolute values. 

Using all the above, it can also be shown that, to 
0{Suj'^) accuracy, the coefficients are given approximately 
by 



b'^'"> - -z7^^'''MVI^2|Tv/|«il , (A6) 



where j^^^^'^ = j^Suj^-^'''') is a real, dimensionless, arbitrar- 
ily small parameter that represents the detuning from the 
crossing point. 

By Eqs. (|A4[) . the coefficients for the second layer are 
given by 



Ini (f,b) 



jif,b) 



(A7) 



These relations show that ^ 1" at the crossing point, 
in which the mode is identically zero. Thus, the gap 
indeed consists only of the crossing point. 
It follows that 



Ni 



di 



f/b I \^f/b\'^ = 2di7^(n2 + ni) [y^TV^Y 



Co 



do = 



^ — ikidi 



^ — ik\ d\ 



1 



"2/ 



"1 ^ « J_ ^^fcirfi 
"2 ■ 



712 



(A4a) 



1 + — ) 6o 

"-2 



(A4b) 

By Eq. (HH), it follows that kidi = tt - ^^5uj = 
k2d2/s = TT - ^^w, where 5uj ^ Sujf/h (sec Eq. p}]) ). 



and similarly. 



Nf/b 



= 2^27^ — ('^2 + ni) (Vn2 T Vm)^ ■ 

«2 



Hence, the normalization constants are 



Nf/b^ I n%,^(j{x)\(j)f/b\'^ = 2ni7^di(ri2 + ni) (^"2 T ' 

(A8) 
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and the mode (jAip in the first layer is given by 

"2 - "1 ik,x 



The diagonal terms reduce to 



4>f/biO < X <di 



"2 T 



-ikix 



with similar expressions for the second layer. Note that 
since for dielectric materials, n2 + rii > ^/n2 + \/n\ 3> 
n2 — n\ > ^/ri2 — ^/ni, the forward (backward) plane wave 
component is dominant for the forward wave envelope / 
(backward wave envelope b) , respectively. Note that the 
size of the coefficient is independent of the frequency, at 
least for non-dispersive materials, so that the dominance 
of the forward/backward component is abruptly switched 
across the zero gap. 

The elements of the uni-directional FB modes (jlOja re 
related to the (standard) FB modes (|A1[) through [3g| 



Using Eq. (|A6|) gives 



from which it follows that 



V'^(a;o) 



1 



^/nphc{xo)(|)ixo) T 



ic 



uj^/nphc{xo) dxo 



(A9) 

It is straight-forward to verify that satisfy a normal- 
ization relation similar to (jAip . 



Appendix B: Calculation of coefficients in envelope 
equations 

In this section we compute the coefficients of the en- 
velope equations (|22ap - (|22cp using the approximate an- 
alytical expressions obtained for the FB modes in Ap- 
pendix [Xj In particular, below we rely on Eqs. (|Aip . 
(|A6l) - (|A7l) and (|A8|) - (|A9|) . All the analytical expressions 
derived below were found to be in perfect agreement with 
numerical evaluation of the definitions (|22aP - (|22cp . 



By Eq. ([^^ and Eq. ^M^, it follows that 



Vab = -—Ifn 



I ^dc^l 1 d<t>, 

-Pfc^ 'Pa^ 

UJa OXo UJb OXo 



dxQ. 



±- 



(Bl) 



Thus, by Eq. (|A5p . we obtain that Vff 



-Vbb = Vg- 



Similarly, one can show that the off-diagonal terms 
vanish, Vab = so that, indeed, there is no reflection- 
coupling in the system. 

9. m(0) 



By Eq. (|22b)) and Eq. ^Mi), it follows that 



'afc 



^{xo)p%(xQ)dxo 

d 



dxo 



p{xo)nphc{xo) (t)*a<t>b 



LOaUJbn%,^(j dxo dxQ 



dxQ. 



Then, it can be shown that 



m 



1°) = pixo)nphcixoWa^bdxo = ^ ^ ^ (B2) 



In Eq. (jB2p . the first term gives the contribution of the zero Fourier component of the modulation and the last two 
terms give the contribution of the ik^Tij Fourier component of the modulation. 

Let us now separate into the diagonal and off-diagonal cases. In the former case, it can be shown that 

™ia = TT^ + TT^-Ti — TT~~7 \ \ f p{xo)cos{2kixo)dx~ I p{xo) cos {2k2Xo) dxA . (B3) 

2ni 2n2 2nidi(n2 + ni) \Jq J^^ I 
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Thus, Eq. (jB3p shows that both diagonal coefficients are identical and that it is only the 2konj Fourier components 
of p{xq) which contribute to the reversal. 

For a piecewisc-uniform modulation, the latter contribution vanishes so that 



(0) (0) "ml , njn2 
TTlfi = mil = . 

'''' 2ni 2n2 



(B4) 



For off-diagonal coefficients, it can be shown that 



.(0) 



712 + ni 2nidi 



For a piecewise-uniform modulation, the contribution of the integrals vanishes so that 



(B5) 



(0) _ n2 - fii / n„a _ n,n2 \ 
~ 2{n2 + ni) V ni 112 J 



(B6) 



3. m 



(1) 



By Eq. 



dxQ. 



This can be shown to equal 



(1) (0) , ^ 



p{xo)n{xo)(pa9bdxo~ / — tt^;^ 

2oJaU}b Jo nphcyxo) 0x0 0x0 



dxo 



Then, using Eqs. (jA6P - (jA7p . it can be shown that 



,(0) _L '^i 



NaNb Jo 



"-2 



(B7) 



It follows that the diagonal terms are identical and that for a piecewise-uniform modulation, ttS^^ 
off-diagonal elements in Ct do not contribute to the dynamics. 



m^"', i.e., the 
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